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Extracting mc{Mc) and /d, ^ from the pseudoscalar sum rules 
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I report here on the (first) direct extraction of the running charm quark mass fhdi') from the D- meson sum 
rules, and on the implications of this result for the estimate of the leptonic decay constants /o(^) . The outputs: 
mc(Mc) = (1.08 ± 0.11) GeV, fo ~ (1.52 ± 0.16)/^, fo, ~ (1.75 ± 0.18)/„ and fa = (1.44 ± 0.07)/^ are in good 
agreement with the existing sum rule results obtained using the pole mass. In particular, the result fo ~ /s 
supports early '87 sum rule results [Q], which indicated a huge 1/m correction to the heavy quark symmetry 
expectation. This talk is based on the paper hep-ph/9712386 and updates the discussions given there. 
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1. Introduction 

One of the most important parameters of the 
standard model is the quark masses. However, 
contrary to the leptons, where the physical mass 
can be identified with the pole of the propagator, 
the quark masses are difficult to define because of 
confinement. Some attempts have been made in 
order to define the quark pole mass within pertur- 
bation theory, where it has been shown to be IR- 
finite 1^ and independent from the choice of the 
regularization and renormalization schemes used 
m . More recently, it has been noticed, in the limit 
of a large number of fiavours, that the resum- 
mation of perturbative series can induce a non- 
perturbative term, which can affect the truncated 
perturbative result, and can, then, limit the accu- 
racy of the pole mass determination Q . However, 
a proper use of such a result needs a resumma- 
tion, at the same level of accuracy (which is not 
often the case) , of the Green function of a given 
process involving the pole mass, where some even- 
tual cancellation between the resummed terms of 
the two series may occur (we plan to come back 
to this point in a future publication). One may 
bypass the previous problems, by working, at a 
given order of perturbative QCD, with the run- 
ning quark masses, which are treated like cou- 
pling constants of the QCD Lagrangian (see e.g. 



1^), and where some non-perturbative-like effect 
is expected to be absent. A lot of efforts has 
been furnished in the literature |^ for extract- 
ing directly from the data the running masses of 
the light and heavy quarks using the SVZ QCD 
spectral sum rules (QSSR) |,|. In this note, 
I shall consider a direct extraction of the run- 
ning charm quark mass using the observed value 
of Md = 1.865 MeV, and study its implication on 
the value of the decay constants , which are 
normalized as = 93.3 MeV where the leptonic 
decay width reads: 



An 



Ml 



(1) 



In this respect, this present work is an improve- 
ment and update of the ones in |lO[§ based on 
the use of the perturbative pole mass, where some 
eventual non-perturbative effects induced by the 
resummation of the QCD series are not taken into 
account. For our discussion, we shall use the av- 
erage value of the experimental widths quoted in 
IHl, from which we can deduce: 
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/d, ~ (1.92 ±0.23)/,. 



(2) 
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while, combining it with the most rehable sum 
rule result for the ratio |p^ : 

/d.//d = 1.15 ±0.04 (3) 

one also obtains: 

/z3 ~ (1.67 ±0.24)/^. (4) 

2. The QCD spectral sum rules 

We shall work with the pseudoscalar two-point 
correlator: 

^^{q^) ^zjd^x (0|TJ,(x)jt(0)|0), (5) 

built from the heavy-light quark current: Jd{x) — 
{nic + md)c{ij5)d, and which has the quantum 
numbers of the D meson. The corresponding 
Laplace transform sum rules are: 



£(r) 



7^(r) 



t< 



dt e — liml>z{t), and 

TT 



dT 



log£(r), 



(6) 



where t< is the hadronic threshold. The latter 
sum rule, or its slight modification, is also useful, 
as it is equal to the resonance mass squared, in 
the simple duality ansatz parametrization of the 
spectral function: 
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ImV-sW^ 2flMl,5(t~Ml) + 

71 

"QCD continuum"e(t - tc),(7) 

where the "QCD continuum comes from the dis- 
continuity of the QCD diagrams, which is ex- 
pected to give a good smearing of the different 
radial excitations |^. The decay constant is 
analogous to = 93.3 MeV; ic is the QCD con- 
tinuum threshold, which is, like the sum rule vari- 
able r, an a priori arbitrary parameter. In this 
paper, we shall impose the tc and r stability crite- 
ria for extracting our optimal results |^ The QCD 
expression of the correlator is known to two-loop 
accuracy (see e.g. Q and the explicit expressions 

^At the optimization scale, its effect is negligible, such 
that a more involved parametrization is not necessary. 
^The corresponding tc value very roughly indicates the 
position of the next radial excitations. 



given in j^]), in terms of the perturbative pole 
mass Afc, and including the non-perturbative con- 
densates of dimensions less than or equal to six 0. 
The sum rule reads: 



C{t) 



dt e"*^ ^^(1 - xf 




e-'''^^\, (8) 



with: 



fix) = - + 2Li2(x) + logxlog(l-x) 

- ^log(l/a;-l)-log(l-a;) 

-I- a;log(l/a; - 1) - (a;/(l - a;))loga;, 

04(04) = -M^{dd) ~ {asG^)/l2TT 



CeiO,) - -^gida^.Y^^a^'d) (9) 

It can be expressed in terms of the running mass 
TOc(i^}_f|, through the perturbative two- loop rela- 
tion 



Me(j^) = me(i^) ^ + ( V) ( ^ + ^ ^ 1 ^ ' ^^^^ 



with: 



P2 

/3i 
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Pi 
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.(11) 



*We shall use the corrected mefficient of the quark-gluon 
mixed condensate given in [ |L1| . This change affects only 
slightly the result. We shall also skip the negligible contri- 
bution from the dimension six four-quark and three-gluon 
condensates. Notice that there is some discrepancy on the 
value of the four-quark coefHcient in the literature. 
^It is clear that, for the non-perturbative terms which are 
known to leading order of perturbation theory, one can 
use either the running or the pole mass. However, we shall 
see that the non-perturbative effects are not important in 
the analysis, such that this distinction does not affect the 
result. 
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where 7™ and /? are respectively the QCD (3- 
function and mass anomalous dimension normal- 
ized as: 

n - i " - i 

which, read for n/ flavours js): 

/3i - 



3. Discussions and results 



/32 = - 
71 = 2 



1 / 2 

1 / f9 

- 5f Uf 

4 I 3 



1 /lOl 



72 



6 



(13) 



As discussed earlier, non-perturbative terms in- 
duced by the resummation of the perturbative se- 
ries can affect this relation However, within 
the approximation at which the spectral func- 
tion is given, the use of this perturbative relation 
should provide the correct expression of the spec- 
tral function in terms of the running quark mass. 
In this way, unlike the analysis in the one 

done in this paper is not affected by the eventual 
existence of such non-perturbative terms related 
to the implicit use of the pole mass in the pre- 
vious sum rule analysis. One should also notice 
that, to the order we are working, the expression 
of the spectral function in terms of the running 
and pole masses differ in the as correction, in- 
duced by the overall leading term appearing 
in Eq. (6). Throughout this paper we shall use 
the values of the parameters |H,|2| given in Table 
1. We have used for the mixed condensate the 
parametrization: 

5(da^.^Grd> = Ml{dd), (14) 



where Mg 



(.8±.l) GeV2 
for four active flavours Itetl: 



A = (325 ± 100) MeV. 



We shall also use. 



(15) 



One can also inspect that the dominant 
non-perturbative contribution is due to the 
dimension- four mc{dd) light quark condensate, 
while the other non-perturbative effects remain a 
small correction at the optimization scale, which 
corresponds to r ~ 0.6 ^ 1 GeV^^ and — 6 ^ 8 
GeV2. 
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Figure 1. Behaviour of fjj versus •ffic{v). The 
horizontal band is the experimental domain of /d. 
The theoretical band is limited by the two curves: 
(a) A=0.425 GeV, v^l GeV w t"^/^, > 7 
GeV2, (Jd)i/3(i GeV)=-238 MeV, M§ = 0.9 
GeV2, {a,G^) = 0.06 GeV^ and (b) A=0.225 
GeV, v K Ma = 1.42 GeV, ^ 6 GeV^, 
(dd) 1/3(1 GeV)=-220 MeV, = 0.6 GeV^, 
(a,G^) = 0.08 GeV2. 



Given the experimental value on Md = 1.865 
GeV, we present our results on from the first 
sum rule for different values of the charm quark 



running mass evaluated at p 



in Fig. 1. 



The second sum rule gives the prediction on Mjj 
for each value of the charm mass (Fig. 2). The 
theoretical band is limited by the two extremal 
values of the QCD parameters used in Table 1. 
Notice that the effect of the errors of the differ- 
ent input parameters is much smaller in the ratio 
of sum rule TI{t). The previous analysis leads to 
the prediction from a two-loop calculation in the 
MS scheme: 



fficiMc) = (1.08 ±0.11) GeV, 



(16) 



where one should notice that, within the present 
accuracy of the data on fu^ , the result comes 
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Table 1 

Different sources of errors in the estimate of fo 



Sources IMfo/Dl 



A = (325 ± 100) MeV 0.12 

(1.20±0.22) GeV 0.08 

mc{Mc) = (1.20 ± 0.05) GeV 0.05 

tc = (6.5 ± 0.5) GeV2 0.04 

(Jd) 1/3(1 GeV) =-(229 ±9) MeV 0.02 

M2 = (0.8±0.1) GeV2 0.02 

(a^G^) = (0.07 ± 0.01) GeV^ 0.01 

Total 0.16 



mainly from the one used to reproduce the ex- 
perimental value of Mo- The radiative correction 
increases the value of ific{Mc) by 10%, whilst the 
non-perturbative term tends to decrease its value 
by a small amount of about 3%, which indicates 
that the running mass is mainly of the perturba- 
tive origin. We have varied the subtraction scale 
L' in the range from t~^^^ « 1 GeV to Mc « 1.42 
GeV, which gives a significant effect (see Table 1) 
due to the log z/ term appearing in Eq. ([l^) . Such 
effect can be eliminated by choosing the subtrac- 
tion point &t v ~ Mc, at which, we shall extract 
the results in this paper. One should also notice, 
in Table 1, that the effect of tc on the result is rel- 
atively small from the value tc — 6 GeV^, where 
one starts to have a r stability until tc> 7 GeV^ , 
where one has tc stability. The main sources of 
errors are due to A and i^. Within the errors, the 
present result is in good agreement (though less 
accurate) with the value |l3| : 

ThciMc) = (1.23iEJ;°^) GeV (17) 

obtained, within the same two-loop approxima- 
tion, from Mj/^. Inversely, we can use the com- 
bined value: 

mc{Mc) = (1.20 ± 0.05) GeV (18) 

from Mjj and from Mj/^ systems on the curve 
fjj as function of fhc given in Fig. 1. Then, one 
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Figure 2. Behaviour of Mo versus ific{v). The 
horizontal band is the experimental value of Mo- 
The theoretical band is limited by the two curves: 

(a) the same as (a) of Fig. 1 but v=l GeV and 

(b) the same as (b) of Fig. 1 but v=lA2 GeV. 



can deduce: 

(1.52 ±0.16)/,, (19) 

where, as can be seen in Table 1, the errors in this 
determination come mainly from the perturba- 
tive parameters: A, mc and the subtraction scale 
V, and, to a lesser extent, from tc and the non- 
perturbative terms. One should notice that the 
radiative as correction and the non-perturbative 
terms have respectively increased the value of fo 
by about 25% and 10%. We can consider that this 
result is an improvement of the previous value: 

/d ^ (1.35±0.07)A, (20) 

obtained by using the perturbative pole mass of 
the charm quark propagator However, the 

good agreement (within the errors) between the 
two results in Eqs ( p^ ) and (|2^) may be an in- 
direct indication that the pole mass defined at a 
given order of perturbation theory, can provide 
a good description of the physical process in this 
channel, and that the eventual non-perturbative 
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power corrections induced by the resummation of 
the perturbative series remain small corrections. 
This observation can also be supported by the 
agreement of the value of the perturbative pole 
mass obtained here and of the one from the J/ tp 
systems, where both values have been obtained 
at two-loop accuracy. A further support of this 
argument can also be provided by the agreement 
of the pole mass deduced from Eq. ( [lO| ) using the 
value of the running mass, with the one extracted 
directly in jlSj . Finally, using the previous value 
of the ratio fnj Id given in Eq. (||), one obtains: 

~ (1.75 ±0.18)/^, (21) 

which is in a better agreement with the data 
quoted in Eq. (H) . A natural improvement of the 
analysis done in this paper is a much more precise 
measurement of fo, or/and an evaluation of the 
QCD two-point correlator to order a^, where, to 
that order, one can have a much better difference 
of the perturbative QCD series by the use of ei- 
ther the running or the pole quark mass. Both ex- 
perimental and theoretical projects are expected 
to be feasible in the near future. 

4. Extension of the analysis to fs 

We have extended the analysis to the case of 
the _B-meson. For this purpose, we use the two- 
loop value |l|,|l|l: 

TOfc(Mfc) = (4.23 ± 0.05) GeV (22) 

from the T systems. We obtain the t and tc sta- 
bihties respectively at 0.2-0.3 GeV"^ and 40-50 
GeV^, similar to the case of the pole mass p0|,p]|. 
Then, we obtain: 

/b~(1.44±0.07)A, (23) 

where the errors come mainly and equally from 
A and mb. The result is in good agreement with 
the analysis where the perturbative pole mass is 
used: 

/b^(1.49±0.08)A. (24) 
This result supports the early '87 sum rule result 


/b « /d, (25) 
indicating a huge \/m correction to the heavy 
quark symmetry expectation fs ^ I/a/Mb. 



5. Extension of the analysis to Jd* 

We have also extended the previous analysis to 
the case of the vector current having the quan- 
tum number of the D* . In this case, and working 
with the correlator having the same dimension as 
V'5 (9^) 07 we do not obtain a t stability as amonge 
other things, the coefficients of the chiral conden- 
sates are opposite of the pseudoscalar ones. Then, 
we conclude that, from this quantity, we cannot 
have a good determination of fo' , within our ap- 
proximation. 
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